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Abstract 

We demonstrate that a suitable coupling between a quintessence 
scalar field and a pressureless cold dark matter (CDM) fluid leads to 
a constant ratio of the energy densities of both components which is 
compatible with an accelerated expansion of the Universe. 
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1 Introduction 



There is a growing consensus among astrophysicists that we live in an accel- 
erating Universe. On the one hand, high-redshift type la supernovae (SNIa) 
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are significantly fainter than expected in a decelerating model (such as the 
Einstein-De Sitter) Q. Although the statistics is still low and extinction 
by interstellar dust may partly account for their low brightness and no con- 
clusive model of evolution of SNIa and their progenitors is still available, 
the acceleration scenario is gaining further ground Q. On the other hand, 
while measurements of the average mass density of the Universe systemat- 
ically fall below the critical density, about 0.3 or 0.4 in critical units (see 
e.g. Q and references therein), the position of the first acoustic peak in 
the temperature anisotropy power spectrum of the CMB strongly suggests 
that the total energy density is critical or near critical [|j. Combining both 
results one may rule out a flat matter-dominated universe (with Q,m = 1 
and S7a = 0) as well as an open universe with no cosmological constant 
(p.M = 0.3 and = 0) at high statistical level ||. More generally, one 
is led to conclude that very likely (i) about two third of the energy of the 
Universe is "dark" (i.e., non-luminous and not subject to direct detection 
via dynamical methods), and (ii) connected to this exotic and elusive energy 
must be a negative pressure, able to violate the strong energy condition. 

The immediate candidate for such exotic energy, a small cosmological 
constant A, poses however an embarrassing question: Why the energy den- 
sity in cold dark matter (which in the absence of interactions redshifts as 
a -3 , where a(t) is the scale factor of the homogeneous and isotropic met- 
ric) and the constant energy associated to A are of the same order precisely 
today? For this to occur one must have fine-tuned initial conditions right 
after the inflationary epoch. This constitutes the so-called "coincidence 
problem" ||. To overcome this hurdle it was suggested that a nearly ho- 
mogeneous but time depending scalar field with negative pressure should 
replace A. This peculiar field, widely known as "quintessence", was inde- 
pendently introduced by Ratra and Peebles |7j and Wetterich ||] well before 
the supernovae results were even suspected. Today a host of quintessence 
models are known both in the realm of general relativity (see e.g., 0], |10[|) 



and in scalar-tensor theories [11]. 

The target of this letter is to clarify a specific aspect of the coincidence 
problem, namely to present an attractor type solution of the two-component 
dynamics which is characterized by a constant ratio of order unity of the 
energy densities of the CDM and quintessence components and at the same 
time admits an accelerated expansion of the Universe. The basic ingredient 
of the corresponding model is to assume a coupling between CDM and the 
quintessence scalar field. It is this assumption of an interacting quintessence 
component by which our analysis differs from most investigations in this field 
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which assume an independent evolution of CDM and the scalar field. "Cou- 
pled quintessence" models have been shown to be useful in handling the 
coincidence problem by Amendola et al. JHJ. While the models of these 
authors assumed a specific coupling from the outset, our strategy here is 
different. We do not specify the coupling from the beginning. We deter- 
mine its structure from the requirement that it shall admit a solution for 
the dynamics of the two-component system of CDM and quintessence with 
a constant ratio for the energy densities. This strategy seems legitimate 
since there does not exist any microphysical hint on the possible nature of 
a coupling between CDM and quintessence. It will provide us with a trans- 
parent phenomenological picture of the "final state" of the cosmic dynamics 
(for a less bleak eschatological scenario see leaving open, of course, 

the question of how this state is approached and whether or not our current 
Universe has already reached it. 



2 Scalar field plus cold dark matter 

Let us consider a two-component system with an energy momentum tensor 

Tik = puiu k + ph ik , (1) 

where h ik = g ik + UiU k and 

P = PS + PM , P = PS + PM ■ (2) 

The subscript S refers to the scalar field component, the subscript M to 
the matter component (i.e. CDM). The energy density and pressure of the 
scalar field are 

PS = \ft + V{<j>) and p s = \ft - V(<p) , (3) 

respectively. The splitting (|2|) implies that there is only one 4-velocity, 

v i _ „< _ „i _ 9 l ^,j 
u ~ U M ~ U S ~ ' 



((f) ^ a is assumed to be timelike.) We postulate that the components do not 
evolve independently but that there exists some interaction between them, 
described by a source (loss) term 5 in the energy balances 

p M + 3i? (pm +Pm) = $ , (4) 
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and 

PS + 3# (ps + PS 
The last equation is equivalent to 

+ 3H<j> + V 



(5) 
(6) 



As already mentioned, we will not specify the interaction from the outset but 
constrain 5 by demanding that the solution to (|J) and (|5|) be compatible with 
a constant ratio between the energy densities pm and ps- It is convenient 
to introduce the quantities I1m and lis by 



S = -3HU M = 3HU S 
with the help of which we can write 

Pm + (pm + Pm + n M ) 

and 



(7) 



(8) 



ps + 3H( Ps + ps + 'ns) = . (9) 

The rewriting of Eqs. (Q) and (||) into Eqs. (||) and (S), respectively, makes 
the dynamic equations formally look as those for two dissipative fluids. The 
fact that there is a coupling between them has been mapped onto the relation 
Hm = —lis between the effective pressures IIm and lis. Some early models 
of power law inflation also share this feature (see e.g., JT^]). Below we shall 
map the interaction term 5 onto a corresponding interaction potential. 

3 Attractor solution and cosmological dynamics 

Consider now the time evolution of the ratio PmI Ps-> 



(10) 
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By introducing the shorthands 
Pm + Pm 



7m 
we obtain 



1 + 



Pm 

Pm 
Ps 
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and 7s 



1M ~ 75 + 
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(12) 
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Obviously, there exists a stationary solution (pm/ps) = for 

n M = -lis = PMPS (75 - 7m) • (13) 
PM + ps 

Since the CDM behaves as dust, i.e. pu <C Pm, we find 



n M * - 1 

or, by virtue of i^ 2 — V = ps ~ p, 



§0 2 + F 



PsPm 
P 



1:2 p 
i^ 2 + V P P 



The coupling term corresponding to this is 

12 



5 = 3H 

or, equivalently, 



1 



U 2 + v. 



PSPM 
P 



(14) 



n M ~ YTo~~Tr~P M = ~P M ■ ( 15 ) 

2 ■ 



(16) 



5 = -3H^p M = -3H ( 7S - 1) -BSE!!-. (17) 

P PS + PA/ 

Introducing the notation r = pm/ Ps = const we may further write 

5 = -3H( ls -l)-^-, or 5 = -3H( ls -l)^—p s . (18) 
r + 1 r + 1 

Invoking the Priedmann equation valid for universes with spatially flat sec- 
tions, 

3H 2 = 8ttG [p s + pm] , (19) 
we have 3H = ^/24nGp, and, consequently, 

<5 = -v / 24^G( 75 -l) / SPM . (20) 

VPS + PM 



With (|lq), in a spatially flat universe equivalent to (|20|), we have identified 
the interaction between the pressureless fluid (CDM) and the scalar field 
(quintessence) that guarantees a constant ratio r of the energy densities of 
both components. 
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To study the stability of this stationary solution against small perturba- 
tions we introduce the ansatz 



PM 

PS 



PM 



PS J st 



+ e 



into (|12|) -the subscript st is for "stationary" . The result is 



m 



3H 
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PS 
PM 
PS 
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PS 
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st 



PS 



PS_( 1+ (PM 



PS 



+ e 



st 



n 



M 



PM 



(21) 



The behavior of the perturbed solution depends on the ratio TLm/pm- For 
the stationary solution itself we may read off 11^ from ([?]) and (|l8|). How- 
ever, for deviations from stationarity an additional assumption is necessary. 
At first sight the most obvious choice seems to be |IIm| oc pu also in the 
vicinity of the stationary solution. As to be seen from (Q), the coupling term 
becomes asymmetric with respect to pu and ps under such conditions. It 
will turn out that a more appropriate choice is the assumption II a/ = —cp, 
where c is a constant c > 0. This type of interaction is symmetric in pu 
and ps- Up to first order in e we find in such a case, 



3Hc- 



-e . 



(22) 



This implies that the stationary solution is stable for r < 1, which is clearly 
compatible with the presently favored observational data pu ~ 0.3 and 
Ps ~ 0.7. Consistency with Um from (0) and (|i~8| ) fixes c: 



1-75 
(1 + r) 2 



(23) 



The positivity of c is guaranteed for 75 < 1. 

With p ~ ps today, the stability condition corresponds to [cf. Eq. (^)] 



P 
P 



n 



M 



PM 



< 



(24) 



Since we seek accelerated expansion, the total pressure p ~ ps must be 
negative, i.e., the potential term must dominate the kinetic term, equivalent 
to 7s < 1. It is remarkable that according to ( p3|) this coincides with the 
condition for c to be positive. From ([/]) and (^) we find that a value 75 < 1 
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implies II m < and 5 > 0. There is a tranfer of energy from the scalar 
field to the matter, which reminds of decaying vacuum energy approaches 
for the dynamics of the early universe (see, e.g., |20|j ). The stationary epoch 
IIm/pm = p/p has to be approached in such a way that 

I XT n-/ 1 Ipl . , 

1 < — . (25) 



Pm P 

Since \TLm\ is proportional to 5, this means, the interaction may be small as 
long as the system is still far from the attractor solution. 

It is expedient to emphasize that the apparently subtle point to assume 
|IIu| oc p instead of \ILm\ oc Pm is essential for the stability properties 
of the stationary solution. Namely, similar considerations as those leading 
to ( |2"2| ) show, that there does not exist a stable solution with accelerated 
expansion for | Hjv^r | oc pu- Therefore, a dependence |IIm| oc p is mandatory 
for a physically sensible solution. This represents a restriction on the type of 
interaction that produces a stationary ratio pm/ps ■ While for the stationary 
solution itself IIm oc pm and TIm oc p are not really different since pM « ft 
the difference becomes crucial if one perturbs the solution. 

Note that the stability is connected to the presence of an effective dissi- 
pative stress in the matter fluid. This parallels the result that the scalar field 
needs the assistence of a dissipative fluid stress for the coincidence problem 
to find solution in spatially flat accelerating Friedmann-Robertsom- Walker 



models 



Given the interaction term (|18|), we may find the dependence of pm and 
ps on the scale factor. Because of pm ~ 0, Eq. (||) with (|l8|) yields 

p M + 3Hp M = -3H ( 7S - 1) , (26) 

r + 1 

while (||) with ( |l8|) results in 

ps + 3H 1S PS = 3H ( 7S - 1) -!— p s . (27) 

r + 1 

Assuming 75, which is in the range < 75 < 2, to be (at least piecewise) 
constant, we obtain 

p s oc a~ v , p M oca"" , v = 3 7S + r . (28) 

r + 1 

Both energy densities happen to redshift at the same rate because we have 
chosen 5 to correspond to the stationary state. With the relationship p oc 
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a v we can solve the Friedmann equation (Jl9|) to find 



- ^T-K)' (29) 

The total energy density redshifts as p oc i -2 , independently of 75 and r. 
Power law accelerated expansion will occur for v < 2, equivalent to 

r + 3 75 < 2 . (30) 

Together with the above derived stability condition r < 1 this amounts to 
75 < 1/3 or ps/ps < —2/3 for accelerated expansion. 
Defining 

8vrG , „ 8vrG . . 

= 3fpP M ' and ^ s = 3^2'° 5 ' ( 31 ) 



we have 



r 1 
= — — 7 , and Q s = — — r , (32) 
r + 1 r + 1 



respectively, and also 
For we find 



Ms = —^Pst 2 • (33) 



1 1 + r 1 
P5 = M ? ' (34) 



Combination with (ITlT) yields 



' 75 (1 + r) 1 1 
6ttG (7s + r) t ' 



(35) 



i.e., evolves logarithmically with time. Furthermore, with the help of (||) 
and (|TT[) it follows that 

Ps = -z = — . (36) 

2-75 75 



which together with (|34|) and (|35| ) leads to 
Since 

V"(0)<£ = F(0) = -2y , (38) 
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by virtue of (35) we obtain 



V(4>) = -WW 



(39) 



where 



and, consequently, 



/ 24vrG 
75 (1 + r) 



(is + r) 



F(^) = Vbexp[-A(^-0 o )] • 



(40) 



(41) 



By similar steps one shows that the interaction term 5 in Eq. (^), given by 
the second expression in (|i~8|), may be mapped onto an interaction potential 
Vint' 

$ - - 2r 1 - 75 



; — Vi n i 



Vi 



■kit 



75 + r 2 - 75 



v{4>) 



Introducing an effective potential 



V eff = V(<f>) + V int 

the equation of motion for the eft field becomes 



+ 3H<j> + Vi 



eff 







(42) 



(43) 



(44) 



It is rather reassuring (although not unexpected, cf.([15|) to find a potential 
( f43|) with (41) and (f42|), substantially backed by some field theories. It ap- 
pears for instance in N = 2 supergravity [|l(|. Likewise, linear combinations 
of exponential potentials naturally arise in theories undergoing dimensional 
compactification to an effective 4-dimensional theory; it is reasonable to 



expect that one of them will eventually dominate [17]. 

With the help of (|30[) the condition for accelerated expansion becomes 



A 2 < 24vrG 



(1-7S) 2 
(l + r)7,S 



(45) 



This is similar but not identical to conditions which have been obtained for 
corresponding solutions in the non-interacting case j|, [l^, 1£] or for different 



types of coupling [10, 13, 18, 21|. These authors started with an exponential 
potential in which A is a free parameter initially. Then they investigated the 
parameter range for which there exists an attractor solution which is also 
inflationary. Our strategy is different insofar, as we have first constructed 
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a solution with the required properties and then read off the corresponding 
parameter combination. 

Notice also that the way the attractor is approached remains open (only 
that in order to guarantee stability the approach, according to (|25|), has to 
proceed from a smaller coupling than given by the stationary solution itself). 



4 Discussion 

We proposed a coupling 5 (given by (|T7|), (|18|), or ([42]) with ([41])) between 
a quintessence scalar field and a CDM fluid that leads to a stable, constant 
ratio for the energy densities of both components, compatible with a power 
law accelerated cosmic expansion. This interacting quintessence approach 
indicates a phenomenological solution of the coincidence problem that af- 
flicts many attempts to cope with late acceleration (especially those based 
in a cosmological constant). Unlike other approaches the potential is not an 
input but derived from the coupling. It remains to be seen to what extent 
this potential is consistent with measurements of the supernovae distances 



[ p2[ once the SNAP satellite comes up with enough SNIa statistics [23|. Al- 
ternative and possibly earlier available tests rely on the Alcock-Paczyhsky 
test for quasar pairs [cf. Ref. 

While focusing on the stationary solution straightforwardly provides us 
with an expression for the interaction which realizes a corresponding state, 
we mention again that this procedure leaves open how this interaction is ex- 
actly "switched on" in order to account for the necessary transition from the 
era of decelerated expansion to that of accelerated expansion. The coupling 
should be ineffective until the condensation of protogalaxies has entered the 
non-linear regime. In a sense, this feature reminds of the "exit problem" 
of many inflationary models. There are attempts to tackle this problem 
with the help of a specific coupling function between (ft and CDM together 



with a separately postulated exponential potential [10|. However, a really 
satisfactory solution is still missing. What one would like to have is an inter- 
action which is negligible in the matter dominated era and asymptotically 
approaches ( |l7| ) for large times. We hope that our stationary solution will 
give an indication for a quintessence-CDM coupling that, aside from char- 
acterizing the stationary state of the late accelerated expansion, smoothly 
joins the previous matter-dominated era of decelerated expansion when one 
goes backward in time. 
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